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Abstract 

We write a 3D equation for three fermions by combining the three two-body potentials obtained by 
the reductions of the corresponding two-fermion Bethe-Salpeter equations to equivalent 3D equations, 
putting the spectator fermion on the mass shell. In this way, the cluster-separated limits are still exact, 
and the Lorentz invariance / cluster separability requirement is automatically satisfied, provided no 
supplementary approximation, like the Born approximation, is made. The use of positive free-energy 
projectors in the chosen reductions of the two-fermion Bethe-Salpeter equations prevents continuum 
dissolution in our 3D three-fermion equation. The potentials are hermitian and depend only slowly 
on the total three-fermion energy. The one high-mass limits are approximately exact. 

In view of a possible perturbation calculation, correcting the remaining discrepancies with the 
three-fermion Bethe-Salpeter equation, we succeeded in deriving our 3D equation from an approxi- 
mation of the three-fermion Bethe-Salpeter equation, in which the three-body kernel is neglected and 
the two-body kernels approached by positive-energy instantaneous expressions, with the spectator 
fermion on the mass shell. The neglected terms are transformed into corrections to the 3D equation. 
A comparison is made with Gross' spectator model. 

PACS ll.10.Qr Rclativistic wave equations. 

PACS 11. 10. St Bound and unstable states; Bethe-Salpeter equations. 
PACS 12.20.Ds Specific calculations and limits of quantum electrodynamics. 

Keywords: Bethe-Salpeter equations. Salpeter's equation. Breit's equation. 
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1 Introduction 

In the treatment of the three-body problem, there is a large gap between nonrelativistic quantum mechan- 
ics (Schrodinger equation) and relativistic quantum field theory (Bethe-Salpeter equation |], || ) . Starting 
from the Schrodinger equation one can of course replace the free part of the hamiltonian by its relativistic 
form. Starting from the Bethe-Salpeter equation, one can try to eliminate the two relative time variables 
to get finally a Schrodinger equation with a compact potential plus a lot of correction terms of various 
origins. Instead of the Schrodinger equation, one can use the Faddeev equations. These equations can 
be derived from the Schrodinger equation, or, in a more general form, from the Bethe-Salpeter equation. 
Faddeev equations for the transition operator give the various scattering and reaction matrix elements 
and the poles (to be computed by basically nonperturbative methods) of this operator in the total energy 
give the spectrum of the three-body bound states. 

Schrodinger's equation is "relatively easy" to solve, but this zero-order approximation does not reflect 
several important properties and symmetries of the studied physical system. These could be recovered, 
in principle, by incorporating higher-order contributions (in general an infinity of them). 

As an intermediate step between nonrelativistic quantum mechanics and quantum field theory we 
shall search for a 3D equation built with the sum of the relativistic free hamiltonians plus two-body 
and perhaps also three-body interaction potentials. Such equations are closely related to the systems of 
coupled Dirac or Klein-Gordon equations of constraint theory [|[ || [| ||, [7), which are still 4D equations 
but exhibit a simplified dependence in the relative times, which could be completely eliminated to get 
a single 3D equation. Besides containing relativistic free hamiltonians, we shall try to make our 3D 
equation satisfy the following list of requirements: 

— Correct nonrelativistic limit. 
- Lorentz invariance and cluster separability. It is always possible to render an equation Lorentz 
invariant by working in the general rest frame (center of mass reference frame) and by building invariants 
with the total 4-momentum vector, although the result may be unelegant and artificial. The Lorentz 
invariance requirement becomes a tool when combined with the cluster separability requirement: when 
all mutual interactions are "switched off", we must get a set of 3 free Dirac equations. This total 
separability can easily be obtained by using as hamiltonian the sum of three free Dirac hamiltonians and 
interaction terms. The real difficulty appears when only the interactions with fermion 3 (for example) 
are switched off. If we want a full cluster separability, the resulting equation for the (12) cluster can not 
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refer to the global center of mass frame anymore, as the momentum of fermion 3 enters in the definition 
of this frame. 

Lorentz invariance and cluster separability can be explicit or implicit (via rearrangements) . The best 
known example of implicit Lorentz invariance is a free Dirac equation solved with respect to the energy: 
it becomes explicitly covariant by multiplication with the P matrix. Other implicit Lorentz invariances 
are not that trivial. For example, the 3D reductions of a Bethe-Salpeter equation are implicitly covariant, 
provided the series generated by this reduction is not truncated. 

— Hermiticity and total energy independence of the interaction terms. In the two-body problem, 
nonhcrmitian interaction terms can be hermitian with respect to a modified scalar product, or made 
hermitian via a rearrangement of the equation. In the three-body problem these rearrangements could 
be more complicated. The hermiticity and the independence in the total energy are linked features, as 
one of these is often achieved at the expense of the other one. Energy dependending interaction terms 
destroy some of the advantages of the use of an hermitian hamiltonian, such as the mutual orthogonality 
of the solutions, and leads to modify the usual perturbation calculation methods. The 3D potentials 
deduced from field theory are generally energy dependent, at least in the higher-order terms. We shall 
require hermiticity and energy-independence (or slow energy dependence) in the lowest order terms at 
least. 

— Correct heavy mass limits. When the mass of one of the fermions becomes infinite, its presence must 
be translated in the equations by a potential (Coulombian in QED) acting on the other fermions. For the 
two-body problem the "one-body limit" in QED is indeed a Dirac-Coulomb equation or a rearrangement 
of it (for example when a projection operator is introduced in order to avoid continuum dissolution - see 
below). In the Dirac-Coulomb equation, the Coulomb potential is already given by the limit of the Born 
term, as the higher-order crossed and ladder terms cancel mutually at the one-body limit. In contrast, 
the one-body limits of the rearrangements contain contributions from all terms || ^, O, [l^] . 

— Perturbative approach. It must be possible to start with a manageable approximate equation and 
indefinitely improve the approximation of the measurable quantities (with respect to the uncalculated 
predictions of the here assumed exact Bethe-Salpeter equation) by adding higher-order contributions. 

— Solution of the continuum dissolution problem. In the relativistic equations for several relativistic 
particles, the physical bound states are degenerate with a continuum of states combining asymptotically 
free particles with opposite energy signs. This often neglected fact forbids the building of normalizable 
solutions in the N > 2-body problem (including the two-body plus potential problem. In the pure two- 
body problem the mixing is prevented by the conservation of the total momentum). The usual solution 
consists in including positive-energy projection operators into the zero-order propagator p3l |l4] , |l5| , 

The modified equations must of course continue to satisfy the other requirements, like the Lorentz 
invariance / cluster separability requirement. 

Section 2 is devoted to the two-fermion problem, revisited in order to define the notations and present 
the building blocks and cluster-separated limits of our future three-fermion equations. In section 3 we 
write directly a cluster separable 3D equation, by combining the relativistic free hamiltonians and the 
3D potentials obtained in the reduction of the two-fermion Bethe-Salpeter equations, putting in each 
3D potential the spectator fermion on the mass shell. Switching off two of the three mutual interaction 
potentials gives then a free Dirac equation for the spectator fermion plus a two-body equation, equiv- 
alent to the two-fermion Bethe-Salpeter equation, for the two interacting fermions. This equivalence 
insures the Lorentz invariance / cluster separability property. Positive-energy projectors solve the con- 
tinuum dissolution problem. Section 4 is devoted to the computation of the heavy-mass limits of the 
three-fermion equation, which are compared with the two-fermion in an external potential equations. In 
section 5, we show that our 3D equation can be obtained from the three-fermion Bethe-Salpeter equation, 
by neglecting the three-body kernel and replacing the two-body kernels by positive-energy instantaneous 
approximations, equivalent at the cluster-separated limits. The neglected terms are thus known at the 
Bethe-Salpeter equation level. A comparison is made with Gross' spectator model jl7| This section 
ends with an attempt of transforming the neglected terms in the Bethe-Salpeter equation into correc- 
tion terms to the 3D equation, in view of a possible perturbation calculation. Section 6 is devoted to 
conclusions. 

2 The two-fermion problem. 
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2.1 Notations. 

We shall write the Bcthe-Salpeter equation for the bound states of two fermions Q as 

$ = G A"$, (1) 

where $ is the Bethe-Salpeter amplitude, function of the positions x\,x 2 or of the momenta pi,p 2 of the 
fermions, according to the representation chosen. The operator K is the Bethe-Salpeter kernel, given 
in a non-local momentum representation by the sum of the irreducible two-fermion Feynman graphs. 
The operator Go is the free propagator, given by the product G01G02 of the two individual fcrmion 
propagators 

_L ^ _ P10 + K 

PiQ -hi + iehi 1 pf - mf + it 

where the h; are the Dirac free hamiltonians 



G °* ~ — u , 77T & ~ J2. 1 ™2 * „•, Pi ( 2 ) 



hi=di.pi + Pirrii (i = l,2). (3) 

The Bethe-Salpeter kernel K should contain charge renormalization and vacuum polarization graphs, 

while the propagators Go; should contain self-energy terms (which can be transferred to K [jl9], |2(J ]). 

In this work, we consider only the free fermion propagators in Go; and the "skeleton" graphs in K, and 
we hope that the inclusion of the various corrections would not change our conclusions. 
We shall define the total (or external, CM, global) and relative (or internal) variables: 

X = -{ Xl +x 2 ), P = Pl + P2 , (4) 

x = X!-x 2 , P= -xipi-p-i)- (5) 
and give a name to the corresponding combinations of the free hamiltonians: 

S = h 1 + h 2 , s = -(h x -h 2 ). (6) 
We know that, at the no-interaction limit, we shall have to get a pair of free Dirac equations: 

(pio - fci)* = 0, (p 20 - = 0, (7) 
where depends on xi,x 2 . Let us also write their iterated version 

(p? - = 0, (p 2 20 - E 2 2 )V = (8) 

with 

E l = ^=(tf + mf)i. (9) 

Interesting combinations can be obtained from the sum and differences of the equations @) or of the 
iterated equations (||): 

(P o -S)* = 0, (po-a)* = 0, (10) 

ff * = 0, (pa - A*)* = (11) 

with 

H„ = 2{(p\ - ml) + (pi - m*)]„,„ = P 2 - 2(E'i + Ef) + V \ (12) 

"-^-^-■m^-^-w (13) 
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2.2 3D reduction of the two-fermion Bethe-Salpeter equation. 

The free propagator Go will be approached by a carefully chosen expression G$, combining a constraint 
like S(po — fi) fixing the relative energy, and a global 3D propagator like — 2itt(Pq — S) -1 ^^- The 
argument of the 8 and the inverse of the propagator should be combinations of the operators used in the 
free equations (at last approximately and for the positive-energy solutions) . Thcr exists an infinity of 
possible combinations §, || 0, 0, §| || || ||, || ||, ||, 0, || |1 @. The best choice 
depends on the quantities one wants to compute (energy of the lowest state, hyperfine splitting, recoil of 
a nucleus, etc..) and on the properties one wants to preserve exactly in the first approximation (cluster 
separability, Lorentz invariance, heavy mass limits, charge conjugation symmetry...). All choices would 
be equivalent if all correction terms could be computed, but this is of course impossible. 

We shall see below that the reduced wave function (from which the relative time-energy degree of 
freedom can be trivially eliminated) is given by 

* = G s Go (14) 

The first choice to be made in the 4D — > 3D reduction is that of the constraint which fixes the 3D 
hypersurface (po = M f° r example) on which we want to work. The remaining of Gs is a purely 3D 
operator, the different choices of which result in different 3D operators applied on a common basic 'J, 
and in different rearrangements of a common reduction series giving the 3D potential. The various 
3D reductions of the literature can thus be classified according to the constraint they use. Once this 
constraint chosen, we can only write different equivalent forms (or sometimes projections) of the same 3D 
equation. It is the unavoidable truncation of the reduction series which makes the difference (numeric, if 
we simply want to compute predictions of field theory; more fundamental, if we want to write constraint 
theory equations with QCD inspired potentials). 

Two natural choices for the constraint are S(po— s), based on the first-order equations (|l0|) and S(po—fJ-), 
based on the second-order equations (|Tl|). In the pure two-fermion problem, this second choice has the 
advantage to fix the relative energy in the rest frame to a simple Po— depending number instead of a 
complicated operator. Other constraints could also be chosen, such as that of Gross Q , which puts one 
particle (normally the heaviest) on mass shell. We made a nonexhaustive review in ref. . 

For a given Po, near Pq = S, the free propagator 

G ° = T77- l -TZ — T TT, 1 , T ■ l ftfo ( 15 ) 

2-P0 + Po - "i + - Po - "2 + »e/i2 

considered as a distribution in po, favors the po ~ s (or the po = /i) region. We can thus approach 
G*o(po) with Gs{po) defined by 

G s (po) = 5(p -s)Gs (16) 



or by 
with 



or 



G s (po) = S(p - f ,)Gs (17) 
dpoG {p ) = -^—^ (3,(32 (18) 
T = o( T i + r 2)> T * = -Tff = tt = sign(/ii). (19) 

r = A++-A-, A«=AiAi Af = ^|^- (20) 

This operator r has a clear meaning in the basis built with the free solutions: it is +1 for hi, > 0, -1 
for hi,h% < and zero when they have opposite signs. It comes from the dependence of the po integral on 
the signs of the iehi. The denominator of Gs must also contain an infinitesimal imaginary part, obtained 
by replacing Po by Po + iePo, as usual. The operator r insures a common sign for poi, po2 and Po, so 
that there is no sign ambiguity in the sum of the imaginary parts. 

There exists of course an infinity of possible choices for Gs . They must however be (at least approx- 
imately) identical on the positive energy mass shell Po = E\ + i?2- An obvious simplification can be 
made by replacing r by 1 (its value for the physical free solutions) or by e(Po) (to get the correct value 
for the corresponding antiparticle states too). The merits of the t or no-r choice are the matter of an 
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old debate. The operator r brings apparently useless complications in the pure two-body case. In the 
two-body plus potential problem, however, the generalization of this operator prevents the "continuum 
dissolution" disease (see below). We shall refer to ( |l8| ) as Salpeter's propagator p^] and to this same 
expression without the operator r as the Breit propagator |}5| . The Breit propagator is thus: 

G B = -2m— (21) 

-TO — O 

For a maximum ease of calculation with the zero-order approximation, one can choose the 3D propagator 
of Schrodinger. Many important physical properties are however lost at this approximation and can 
only be recovered by the inclusion of higher-order contributions. By contrast, the 3D propagator of 
Sazdjian, based on the second-order equations (hlh , is more complicated but leads to a covariant Gs 
when combined with the 5(po~- /i) constraint || |12|, jl9| . In the remaining of this section we shall work 
with the approached propagator 

G s = S(po-s)AG B (22) 

where A can be r (Salpeter), 1 (Breit) or A ++ . 

We shall write the free propagator as the sum of the approached propagator, plus a rest: 

G = Gg + Gr. (23) 

The Bethe-Salpeter equation becomes then the inhomogeneous equation 

$ = G K<i> = (G s + G R )K<S> = * + GrK®, (24) 

with 

* = G S K<P (= GsGq 1 ®). (25) 
Solving (formally) the inhomogeneous equation (^J) and putting the result into (|25|), we get 

* = G S K(1 - GrK)- 1 ^ = GsK T f (26) 

where 

K T = K{1 - GrK)- 1 = K + KGrK + ... = (1 - KGrY 1 K (27) 

obeys 

K T = K + KGrKt = K + KtGrK. (28) 
The reduction series (^7|) re-introduces in fact the reducible graphs into the Bethe-Salpeter kernel, but 



with Go replaced by Gr. Equation (|26j) is a 3D equivalent of the Bethe-Salpeter equation. It depends 
on the choice of Gs- 

The relative energy dependence of eq. (|2^) can be easily eliminated: 

* = 5(j> -s)i> (29) 

and ip obeys: 

(30) 

-To — 5 

where V is proportional to Kt with the initial and final relative energies fixed to s : 

V =-2iTT^ 2 K T {s,s). (31) 
In less compact but more precise notations: 

l3i/3 2 K T (s,s) = J dp' dp Q S(p -s)Pi(3 2 KT(Po,Po)S(po-s). (32) 

Note that we write (p' ,po) but (s, s), as we keep s in operator form. 

When A — 1 the operator AV is hermitian. When A = A ++ the operator AV is hermitian in the 
A ++ = 1 subspace, and we can write 

^ =— L-A++VA++V. (33) 
"o — " 



G 



When A = t the operator AV is hermitian in the r 2 = 1 subspace, and we can write 



^=-J_tVt^. (34) 
r — <-> 

Let us recall that the potential V depends on Gg and is thus not the same for the different choices of 
A. 

These 3D reductions can also be described in terms of transition operators. The 4D transition operator 

is 

T = K + K G K + ■■■ (35) 
and the 3D transition operator is, when A = A ++ , for example 

T 3D = A ++ VA ++ + A++VA++ l - A++VA++ + • • • (36) 

P - b + ie 

The relation between these two transition operators is simply 

T 3D = -2in0i02 < T >, (37) 

with 

< T > ee X 2 A++ X 2 T(s, s) A++. (38) 

We have indeed 

^ 0102 T 3D = <K T > + <K T > 5 ~ 2 '^ . 0102 <K T > + ■■■ 

ZlTT Pq — b + 16 

= <K T > + <K T G S K T > + ■■■ = < K T (1 - G s K T y 1 > 

= <K{1-G R K)- 1 {1- GsKil-GnK)- 1 )- 1 > 

= < K (1 - G R K- GsK)- 1 > = < K{1 - GqK)- 1 > = < T > . (39) 

A Bethe-Salpeter equation leading directly to the same 3D reduction can be obtained by replacing the 
kernel K by the instantaneous kernel < Kt > ■ The corresponding transition operator becomes in this 
case 

T = <K t > + <K t >Gq<K t >+---=<T>. (40) 

In contrast with the textbook cases, the potentials deduced from field theory are in general energy- 
dependent (although one starts often with an energy- independent approximation). The eigenvalue of the 
hamiltonian is then a function of the energy itself, and the energy spectrum is given by the solutions of 
the algebric equations 

E = E x (E) (41) 

but the corresponding wave functions arc no more orthogonal with the the usual scalar product, so that 
the usual methods of perturbation calculation have to be revisited |3(| 

Until now we did not specify a reference frame, and our equations are not manifestly covariant under 
the Lorentz group. The consequences of this must however be carefully discussed. 

- Trough not manifestly covariant, equations like ( p6|) can always be transformed back into the original 
Bethe-Salpeter equation, provided the reduction series ( ^7|) is not truncated. The inclusion of higher- 
order terms will thus, in principle, improve an approached covariance. Beyond this somewhat trivial 
consideration, it would of course be interesting to be able to simplify the equations in a way which 
preserves the covariance. If we use the covariant two-body Sazdjian propagator J(| based on the 
second-order equations (|ll|), instead of Breit's, it becomes possible to truncate the reduction series in 
a covariant way (one might for example keep only the first term of the reduction series at the ladder 
approximation, i.e. the Born term). 

- Even truncated, our equations could be made formally invariant by introducing a unit vector n — 
(1,0) (we can call this vector the laboratory time unit) and by making all elements in our equations 
invariant by using scalar products with n. If we want more than a formal invariance, this unit vector can 
however not be external to the system: one must therefore define n as Pj \J P 2 (the time unit of the 
center of mass reference frame), assuming P 2 > 0. 
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- It must be noted that our non-covariant equations admit, among others, P = solutions which are 
identical to the P 2 > solutions of the covariant equations in the center of mass reference frame. We 
could speak of " weak covariance" . 

- An equation like ( |33| ) can be assumed to have been written in the center of mass reference frame 
and made explicitly covariant by using the vector P/V P 2 . This covariant form could be used to write 
the system of equations in other reference frames. There is no reason to do that in the pure two-body 
case, but it will become necessary when other objects (as an external potential or a third particle) are 
present. The Gb 5(po ~ s ) combination written in the center of mass reference frame looses its simplicity 
when covariantized. 

The cluster separability property is clearly satisfied, as the "switching off' of the mutual interaction 
leads to a pair of free Dirac equations. 



3 The three-fermion problem. 

3.1 A cluster separable 3D equation. 

A three-fermion 3D equation, inspired by the two-fermion equation ( |30| ) with A = 1 could be 

i> =77^(^2 + V 23 + V 31 )V (42) 

with S = hi + h-2 + /13. The potentials are the two-body potentials defined by (pl|). Each two-body 
potential Vij was a function of the partial energy PyOj which we shall fix to its (ij)+k cluster-separated 
limit Po — hk. At the V23 = V31 = limit, for example, we shall thus get two completely independent 
equations: 



■Pi 20 — hi — hi 



1pl2 = V 12 (Pl2o)lpl2, P30^3 = h 3 i/j 3 (43) 



with 

1p = ^12 03, Po = P120 + P30- (44) 



Our 3D equation (42) satisfies thus clearly the cluster separability requirement. Furthermore, the three 
cluster-separated limits are exact equivalents of the corresponding two-fermion Bethe-Salpeter equations. 

In the two-body problem all quantities can be defined in the center of mass reference frame. In the 
two-body plus potential problem we had to start in the laboratory frame but to consider also the center 
of mass reference frame at the no-external potential limit. In the three-body problem we must start in 
the three-body center of mass frame (unless we are satisfied with a "weak covariance") and consider the 
center of mass reference frames of the three possible two-body subsystems obtained by cluster separation. 
The fact that the three cluster-separated limits are exact insures in principle the Lorentz invariance 
/ cluster separability requirement: at the cluster separated limits the two-body equation can indeed be 
transformed back into a covariant two-body Bethe-Salpeter equation. There is no necessity of introducing 
Lorentz boosts by hand. 

Our two-body potentials are the sum of an infinity of contributions symbolized by Feynman graphs. 
Keeping only the first one (Born approximation) or a finite number of them renders the Lorentz covariance 
of the two-fermion clusters only approximate. A Born approximation preserving the Lorentz invariance / 
cluster separability property can be obtained by using another 3D reduction based on a covariant second- 
order two-body propagator of Sazdjian |q, |l2, 19 (combined with a covariant substitute of A ++ to prevent 
continuum dissolution - see below). This leads to a 3D three-cluster equation which is covariantly Born 
approximable, but more complicated 

3.2 The continuum dissolution problem. 

Unfortunately, our 3D equation (^) suffers of continuum dissolution. It is indeed possible to build a 
continuum of solutions with any a priori given total energy by combining asymptotically free fermions 
with opposite energy signs. Any physical bound state is thus degenerate with such a continuum and the 
building of normalisable bound state wave functions becomes impossible [|l3[ 14, [lf| [l(|. In the pure 
two-body case the energy of a system (+, — ) in the total rest frame is 



E 1 -E 2 = yjp* + m\ - VF 2 + ml = ' j= 2 (45) 
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and lies thus between mi — 7712 (whichever the sign) and zero. We have thus no problem if we make the 
assumption that the energies of the bound states lie between | mi — m-i | and m\ + m 2 . Below | m\ — 7772 \ 
we meet the "strong field" problem. The strong field problem for the one-body plus potential and two- 
body sustems, the continuum dissolution problem for the two-body plus potential and three-body systems 
are both consequences of the possibility of pair creation. 

In the three-body case, the mixing with asymptotically separated (12) (3) subsystems of opposite 
energy signs is excluded by the total momentum conservation, as in the two-fermion case. For the mixing 
with three-fermion asymptotically free states, let us consider for example 

E 3 -E l ~E 2 = \]pl + ml- yjfi+ml- ^Jp^+mJ, pi + p 2 + p 3 = 0. (46) 

There is no lowest value (we can for example have p 3 = and p\ — —p>2 arbitrarily large) . The highest 
values are obtained when p\ and p% have the same direction, and for pi/mi =77*2/7772 = — 373/(7771+7712). 
In this case we have 

E 3 -Et-E 2 = Jf 3 + m§ - ^p| + (mi+m 2 ) 2 (47) 

so that E 3 — Ei — E2 lies finally between —00 and (771,3 —Tn\ — 7712)^(7713 — mi —1112). Symmetrically 
E\ + E2 — E 3 lies between (mi+m 2 — m 3 )0(rn 3 — rax— TO2) and 00. If we assume that the energies of the 
bound states lie between (7711+7772 + 7773— 2Inf(mj)) (i.e. above the highest negative-energy threshold) 
and (mi+m 2 +m 3 ), there is no degenerescence with the "one plus-two minus" states. On the contrary, 
no weak field assumption could prevent the degenerescence with the "two plus-one minus" states. 

This continuum dissolution problem can be cured by simply introducing the products of noncovari- 
ant projectors A+ + = A+A+ (on the eigenstates of the free hi,hj with positive eigenvalues) into the 
corresponding Gsij, defining the two-body potentials, and around these potentials. This changes finally 
equation ( ^2| ) into 

4, = — L- ( A++V 12 A+ + + A2^2 3 A+ + + A3++F31A++ ) ^ (48) 
Fq — o 

built with the three two-fermion equations (|33|). There exist solutions of eq.(^8|) which are entierely 
contained in the (H — | — h) subspacc. These solutions obey an equation with a spinless Salpeter free term: 

(P -E 1 ~E 2 -E 3 )4,= A+++ ( Vu + V 23 + V31 ) A+++ (49) 

It is also possible to build a three-fermion equation with three Salpeter equations: 

i> = 5 1 g ( T 12 VuT? 2 + T23V A 23r 2 2 3 + T 31 V 31 T^ ) 1p. (50) 
Fq — & 

This equation admits solutions lying entierely in the (+ + +) ® (H ) © ( — I — ) © ( h) subspace. 

The remarks above about the cluster separability and Lorentz invariance of equation ( [42] ) remain true 
for equations ( p8^ and (|50|). 

In the remaining of this work we shall adopt the positive-energy equation ( f48| ) as our basic three-cluster 
3D equation. 



4 Heavy mass limits and external potentials. 

4.1 Two fermions in an external potential. 

The two-fermion plus potential problem can be approached in two ways which we shall call the two- 
body and the three-body approaches. In the two-body approach, the external potential is included in the 
definition of the creation and annihilation operators of field theory. Practically, we can keep the equations 
obtained in the treatment of the pure two-fermion problem, using simply a generalized definition of the 

hi = di .pi + fami + Vi(xi) (51) 

where Vt is the external potential acting on the fermion i. All quantities can be expanded on the basis 
built with the eigenstates of hi and ft 2 . In the three-body approach, the free creation and annihilation 
operators are used and the external potential is treated as an heavy third body. This approach will be 
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presented as an heavy mass limit of the three fermion problem in subsection 4.3. Here we shall adopt the 
two-body approach. 

The equations for two fermions in an external potential must be written in the laboratory reference 
frame (in which the external potential is defined) and are not Lorentz invariant. If we "switch off" the 
mutual interaction we get a pair of uncoupled equations for two independent fermions in an external po- 
tential. If we switch off the external potential the equations remain written in the laboratory frame, which 
still refers to the vanished external potential. If the reduction series is not truncated, this no-external 
interaction limit of the equation is equivalent to a pure two-fermion covariant Bethe-Salpeter equation. 
If the reduction series is truncated, the pure two-body cluster equation is only "weakly covariant". If we 
want a truly cluster separable system, we must then use equations which are covariant at the vanishing 
external potential limit, such as Sazdjian's equations with the generalization ( |5l| ) of the hi || [l^, ^9|. In 
the two-body plus potential problem the total spatial momentum is no more conserved and the mixing 
with the continuum forbids the building of normalizable solutions. A easy way of seeing how this hap- 
pens is to try to build these states by perturbations of the no-mutual interaction equations: the operators 
hi, hi are then diagonal and the higher-order contributions contain denominators in Pq — h\ — h 2l which 
can vanish on a continuum. In Salpeter's equation the operator r kills the matrix elements of r V 
between the r 2 = 1 states (+,+),(—,—) and the r 2 = states (+,—),(—,+) in the basis built with 
the eigenstates of (hi, hi), forbidding the mixing. 



4.2 Heavy mass limits in the two-fermion problem. 

When the mass of one of the fermions goes to infinity, we must find the equation of the other fermion in 
a potential (a Dirac equation with a Coulomb potential in QED). This limit can be obtained directly, or 
via a rearrangement of the equation. With Breit's equation with the second-order constraint and with 
Sazdjian's equation, the higher-order ladder and crossed terms cancel mutually at the one-body limit, so 
that the correct limit of the potential is already contained in the Born term. With Salpeter's equation 
with a second-order constraint, the limit is also a Dirac-Coulomb equation, but solved with respect to 
the A + ip part of the wave function: the A~ip part is transformed into higher-order contributions to the 
potential. The physical content remains identical to that of a Dirac-Coulomb equation, but the correct 
limit of the potential is no more entierely contained in the Born term, so that a truncation of the potential 
would spoil the one-body limit g |, |(J || ||. 

4.3 Heavy mass limits in the three-fermion problem. 

When m3, for example, goes to infinity (two-body limit), we get, writing Pq = Wn + UI3 in our basic 
three-cluster equation: 

W 12 if) = [hi + h 2 + A++ Vu A++ + A+ V+ A+ + A+ V+ A+) V (52) 

with (Pq — h%) replaced by Wn in Vu- The potential V 2 + is given by the series 

V+ = V 2 +V 2 A f V 2 + ... (53) 

Wn - hi - hi 

where V 2 is an external potential acting on fermion 2 (it is a Coulomb potential in QED, if we use the 
second-order constraint (|l7|), but in general it could still depend on (Wi 2 — hi) ). The potential is 
given by a similar formula. 

The differences with the two-body approach of the two-body plus potential problem are: 

- The projectors A^ are the free ones. 

- We have now projectors around the external potential terms V 2 + , . 

- These external potential themselves are now given by the series ([33]). 

- V 2 could still depend on (Wn — hi) and vice- versa. 

The cluster separability property in the three-body way survives the high-mass limit: switching off V 2 
and Vi leads to the equation for the (12) two-fermion system, switching off V 12 an d Vi leads to a free 
Dirac equation for fermion 1 with the A^ projection of the Dirac equation for fermion 2 in the external 
potential V 2 . In the two-body approach of the two-fermion plus potential problem we had however more: 
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switching off only the mutual interaction led to a pair of independent equations for each fermion in the 
external potential. Here, in the three-body approach, the equation does not split perfectly into two 
parts: we can write W12 =W\ + W2, but we have W12 — hi instead of W2 in the series defining V 2 + 
and vice-versa. This is a consequence of the energy dependence introduced into (^3|) by the use of the 
anti-continuum dissolution projectors, combined with our choice of putting the spectator fermion on the 
mass shell in each two-body interaction, and neglecting the three-body terms which could balance this 
modification. The discrepancy is of order V A . At the two-fermion plus potential level, a better but 
still not perfect equation would be obtained by replacing h\ by hi + V\ in the series defining V 2 + and 
vice- versa. 



5 Towards a perturbation calculation. 

5.1 Bethe-Salpeter equivalent of our basic three-cluster equation. 

The three-fermion Bethe-Salpeter equation can be written 

3? = [G01G02K12 + G02G03K23 + G03G01K31 + G01G02G03K123] <& (54) 

where K123 is given by the sum of the purely three-body irreducible contributions. We would like to get 
our basic three-cluster 3D equation by approximating the three-fermion Bethe-Salpeter equation, and to 
possibly recover the neglected contributions in a perturbation calculation afterwards. A three-fermion 
Bethe-Salpeter equation with instantaneous (i.e. independent of the relative energies) and positive-energy 
kernels can be transformed into a 3D Salpeter equation Q . We shall thus neglect the three-body kernel 
and replace the two-body kernels by instantaneous positive-energy kernels (with the spectator fermions 
on their mass shell) which are equivalent at the cluster-separated limits: 

K X 2 « < K T12 {P -ha)>= 0i0 2 A^ 2 + ^12(^0-^3) A++ , • • • . (55) 

ZlTT 

The Bethe-Salpeter equation becomes 

$ = ^ GoiG 02 G 3 01010a [ V12 A++ 12 + A++ ^23 A++ ^23 

+ A++F31A++V31] (56) 

where 

ipij (Pko) = 0k G^ J dpijo $. (57) 
This leads to a set of three coupled integral equations: 

^12(^30) = 7— / dp 12() G01G02 0102 [ Af 2 + V12 A++ feOso) 
2nr J 

+ A++ V 23 A++ V 23 (Pio) + A++ V31 A++ ^ 3 i(P2o) ] (58) 
where pio,P20 must be written in terms of Po,P30,Pi20 '■ 

Pq — P30 , Pq — P30 ,rn\ 

Pw = ^ h Pl20 > P2a = 2 Pl20 ' ^ ' 

and similarly for 023 and "031. We shall now search for solutions Vij^fco) analytical in the lower 
Im(pfeo) <0 half plane and perform the integration ( ^8|) by closing the integration paths around these 
regions. The only singularities will then be the poles of the free propagators. The result is 

A++ 

<MP30)= 7p F — - [A+ + ^ 1 2A++0 1 2(P3O) 

(P - S) - (P3Q - h 3 ) + ie 
+ A++ V 23 A++ ip 23 (hi) + A++ V31 A++ 03i (h 2 ) },■■■■ (60) 
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Using these equations to compute the projections — 2in tpiji^k) we find that these three expressions 
are equal (let us call them ip) and obey the A +++ projection of our three-cluster equation ([48|): 

^ = 4^ [^12 + ^3 + V 31 ] A+++^. (61) 

Furthermore, we have also 

J dp § = J dp 10 dp 2 odp3o S(pio + P20 + Pso - Po) $ = J dpko G k (3k tpij (pko) 

= -2i7rA+Vy(/i fe ) =i/j. (62) 
Solving ( J60| ) with respect to ipi% (P30) gives 

^12(P30) - ^ — ,c A l + !++T/ A+ + ^- [y23 + y3l] A+++1P (63) 

2^7^ (P - P30 ) - (Su + Ai 2 + V12AJ7) + it 
which confirms the above assumption about the analyticity of the 4>ij(Pko)- 



5.2 Relation with Faddeev formalism and Gross' spectator model. 

The nonrelativistic Faddeev equations can be obtained by transforming Schrodinger's equation (in a first 
step, without three-body terms) into a set of three coupled equations for three parts of the transition 
operator ( denotes the contribution of all graphs beginning by a (ij) interaction). The input is the set of 
the three two-body transition operators and the resulting series expansion contains only connected graphs 
(never twice the same two-body transition operator). This formalim is well adapted to the description 
of the various scattering processes, such as (12) + 3 — > (12) + 3 (elastic scattering), (12) + 3 — > (12)* + 3 
(excitation), (12) + 3^1 + (23) (rearrangement), (12) + 3^ 1 + 2 + 3 (breakup). The (123) bound 
states correspond to the poles of this transition operator. 

The structure of these equations can be generalized to relativistic equations which arc not necessarily 
(exactly) reducible to a single 3D equation. 

In Gross' spectator model 0, Faddeev's type equations are deduced from the Bethe-Salpeter 
equation and the transition operator is written in terms of the two-body transition operators. The relative 
time variables are then eliminated by putting, in each three-body propagator, all the " offmassshallness" on 
the only fermion which interacts before and after. The Lorentz invariance-cluster separability requirement 
is satisfied by applying suitable Lorentz boosts on the two-body transition operators. 

The main difference between our approach and Gross' approach comes from the fact that we are 
(presently) interested by the (123) bound states only, to be computed principally by using a single 3D 
equation (such as the one from which Faddeev starts). Instead of working with the two-body transition 
operators, we work with the two-body potentials. We can however present our approach in terms of the 
two-body transition operators. Our approximation of each two-body transition operator is then unique 
and does not depend on the operators which come in front and behind it in the expansion of the Faddeev 
equations. This enables us to describe our model by a single 3D potential equation, by a set of three 
Faddeev equations, or by the expansion of the three-body transition operator in terms of the two-body 
ones. Another option of our model is to not introduce Lorentz boosts by hand: the Lorentz-invariance / 
cluster separability requirement is exactly satisfied if we do not truncate our potentials (practically, this 
means that the satisfaction of this requirement and the approximation of the potentials can be improved 
together). Let us go back to the three- fermion Bethe-Salpeter equation and neglect for simplicity the 
three-body irreducible kernel -^123- Defining 

$12 = G01G02 if 12 • • ■ $ = $12 + $23 + $31 , (64) 

we get 



Writing then 



(1 - GoiG 2^1 2 )$ = $23 + $31 
$ = (1 + G iGo 2 Ti 2 ) ($23 + $31) 
$12 = G01G02T12 ($23 + $31 )■ 

$12 = G01G02G03 Xl2i •'• 



(65) 
(66) 
(67) 

(68) 
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in order to factor out the propagators, we get 

Xl2 = 0102 T\2 G01G02 0102 ( X23 + X31 )• 



(69) 



The (12) transition matrix element corresponding to our approximation (^5|) of K\ 2 is 

T 12(Pl20> P120, Pa -P30) ~ ?1 2(P30) 
= < ^T120Po - M > (1 - G01G02 < Kt1 2 (Po - M >r X - (70) 

This r^feo) is analytical in the Im(p3o) < half plane and 

T° 2 {h 3 ) = <T l2 {P Q -h 3 ) >. (71) 

This approximation, combined with equation (|69|), implies that Xi2(Pi20iP3o) is independent of p\ 2 a (let 
us write thus Xi2(P3o))- Equation (69) becomes then 

Xi2(P3o) = 0i02T° 2 (p 3O ) / dp 12O G Q1 G O2 1 2 [x23(Pio) + X31CP20)] (72) 



where Pw,P2o must be written in terms of Po,P3o,Pi20- Equation (|72|), together with similar equations 
for X23(pio) and X3i(.P2o), admits solutions which are analytical in the Im(pfeo) < half planes. At 
Pko = hk, we get the 3D Faddeev equations 

Xi2(h 3 ) = T? 2 D (P - h 3 ) [ Xaafai) + X3 i(h 2 )}, ■ ■ ■ (73) 

which can easily be transformed back into our basic three-cluster potential equation (by performing in 
the reverse order the transformations made above at the 4D level). Note that we would get the same 
result by approaching directly T\ 2 by T° 2 (h 3 ) instead of Ti 2 (p 3 o). This means approaching if 12 by a 
given function Ki 2 (p 3 o), analytical in the Im(p3o) < half plane and equal to < Kt\2 (Pq ~ h 3 ) > at 
P3Q = h 3 . 

The key point of the manipulations above is the dominance of the positive-energy poles of G01G02 
in ([39|). This was obtained by approaching Ktyi (or T\ 2 ) by a constant with A+ 2 + positive-energy 
projectors. We could try to insure the dominance of these two poles more economically. Let us come 
back to equation ( |69| ) without making any approximation on T\ 2 . The elements of ([39]) are then the 
operator and functions 

T i2 (P120. P120, Po ~ P30) (74) 

X12 (p'i20> P30 ), X23 (P230, PlQ ), X31 (P310, P20) (75) 



and ( p9[ ) must be integrated with respect to puo, with p 3 o fixed. We must thus write P230: P10, P310, P20 
in terms of P120, P30- Searching as above for solutions Xij ( Pijo 1 Pko ) which are analytical in the 
Im(pko) < half planes, we shall close our integration path clockwise (counterclockwise) in front of 
X23 ( X31 ) an d keep thus the pole of G01 ( G02 ), which puts fermion 1 (2) on its positive-energy mass 
shell. The elements of ( |69| ) are then replaced by 

T i2 (p'rao! P120, Po ~P3q) -> T 12 (p' 120 , s i2 T P ° 2 ± P3 ° 2 ^ , -Po - Pao) (76) 



/ 



— 2?7T 

d Pl 2 Q G 01 G 02 M 2 - — — — — — (77) 

(P -S)- (p 3 o - h 3 ) + ie 

P — S 

X23 (P230, PlQ ) -> X23 ( S 2 3 H 2 ( P30 ~~ ^ 3 )' hl ) ^ 78 ) 

Pn — S 

X31 (P310, P20 ) ->• X31 ( S31 2 h ( P30 ~~ ^ 3 )' h2 ) ( 79 ) 



where we take the upper signs in ( |76[) in front of X23, the lower signs in front of X31- 

The manipulations above are submitted to some restrictions on the . In the pko variable the Ty 
must be analytical in the lm(pko) < half planes. In the p'^ and p^o variables the must be 
analytical in the whole complex plane. Moreover, the must also be asymptotically bounded in the 
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three variables and contain Ay positive-energy projectors. These conditions are of course not satisfied 
by the exact transition matrix elements, but we shall assume that the singularities and the negative- 
energy parts of the Tij can be neglected in the computation of the integrals with respect to the relative 
times. 

If we take then the equation for X12 a t P3o — h 3 and p' o = sf 2 = s i2 i 5(^0 — S), etc., we get a 
closed system of six 3D equations: 

X + - -2zn fafa ( T+- x+ + T++ X 3i ) (80) 

Xl2 = -2i7r fcfo ( Tf 2 - ^ZTs + T i2 + X31 ) (81) 

and similarly for X23 an d X3i- 

In each T^ two fermions are on the mass shell: the spectator fermion k and the fermion which is not 
going to interact next left (p^ ) or next right (pijo ). This is of course the philosophy of Gross' spectator 
model 

If we replace s ± by s in the TJy we get our basic three-cluster model. This supplementary approx- 
imation can be introduced by noticing that we already neglected the contributions of the singularities 
of in the p'^ and p^o complex planes. The dependence on these variables can thus not be very 
strong, as we know that a function which is analytical and bounded on the whole complex sphere must 
necessarily be a constant. Using the same argument of the consistency of the approximations, we could 
also argue that Ty constant implies Kxij and constant with Kxij — Kij . Our basic three-cluster 
equation is thus not a priori a worse approximation of the three-fermion Bethe-Salpeter equation than 
Gross' equations, nor a better approximation than the simpler positive-energy instantaneous approxi- 
mation or the corresponding Born approximation. It reflects the choice of preserving exactly the exact 
cluster-separated limits. 

It is thus important to evaluate the importance of the neglected contributions in a perturbation 
calculation. 



5.3 Correction terms. 

The contributions neglected above in the transformation of the Bethe-Salpeter equation into a three- 
cluster 3D equation are given at the Bethe-Salpeter level. We would like to transform them corrections 
to the 3D equation. However, the exact Bethe-Salpeter equation has no more the simple analyticity 
structure in the relative energies which allows the 3D reduction at the positive-energies instantaneous 
approximation (|55|). 

We made several attempts to build a 3D perturbation calculation starting with our basic three-cluster 
equation. Some of them led to interesting results, but not to the kind of 3D perturbation terms we 
wanted. In the two-fermion problem, a perturbation expansion can be built around an approximation of 
the propagator (as in section 2), of the kernel (as in [p8| for example), or of both. Our approximation 
( [55| ) for the three-fermion Bethe-Salpeter equation is an approximation of the kernels: the three-fermion 
kernel is neglected and the two-fermion kernels are replaced by series based on approximations of the 



two-fermion propagators. In order to transform the correction terms, we shall write, instead of (55) 

& (p'30 -P3o) K12 Oi 20 , P120 , Pq -P3o) 

« 6(p' 30 -h 3 )A+ < K T12 (P Q -h 3 ) >= K° 2 (p' 30 ). (82) 

This gives directly the (+++) projection of our basic three-cluster equation for the integral (^32|) of 
the Bethe-Salpeter amplitude with respect to the relative times. The unapproximated Bethe-Salpeter 
equation takes the form 

* = [G iG a2 K° 12 + G 02 Go3^ 2 3 + GosGoiK^ + R] $ (83) 
and leads to the 3D equation 

^ = Jd Po A+++(l - R)- 1 [GoiGo 2 K° 12 + G 02 Go3^ 2 °3 + GqsGoiK^] i> 
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A+++ 
Po-S 



(V 12 + V 23 + V 31 )A +++ + R 



1> 



with 



R = Jd Po A+++ T -^L [G 01 G 02 K° 12 + G 02 G 03 K° 3 + G^G^K^] 

In terms of x = (l — R)"^ '■ 

{Pa-S) X = A+++(V 12 + V 23 + ^i)A+++(l-i?)- 1 x . 



(84) 
(85) 

(86) 



In this last form the correction terms are more symmetric and less total energy dependent. At the 
approximation of positive-energy instantaneous two-body kernels, as ( j55| ) , we expect R (or at least Rip ) 
to vanish. This is necessary if we want it to be of higher-order in the general case, but it is not automatic: 
although the "approximation" ( |8^ ) led to our basic three-cluster equation, equation ( |86] ) could turn out 
to be a complicated equivalent of it (we met this situation in previous promising attempts of perturbation 
calculation). For positive-energy instantaneous two-body kernels, we have 



R — G01G02K12 + Gq 2 Gq 3 K 23 + Gq 3 GoiK 31 = i?i2 + i?23 + R31 (87) 

tf£(P30.P3o) = K 12 [S(p' 30 -p 30 ) - A+S(p' 30 -h 3 )] . (88) 

This implies immediately that the terms ending like Ri 2 GoiGq 2 Ki 2 are zero. For terms ending like 
Ri2Gq 2 Gq 3 K 23 we have to examine the action of four possible operators at left: 



dp (. 



R 



23- 



R3I, Rl2- 



The three first ones lead to contour integrals which _give zero. With R\ 2 we must perform the p\ 2 Q integral 
and then apply again one of the four operators (|39|) at left, etc... This recurrence ends necessarily with 
one of the first three operators ( |S9l ) which gives zero again. 

In the general case our correction terms are thus really higher-order terms. Two aspects remain 
however to be investigated: 

- It seems that our method brings corrections to the two-fermion potentials as well. At the cluster- 
separated limits, our two-fermions equations would then be replaced by equivalent equations, unless we 
could show that these two-fermion correction terms vanish. 

- Our three-fermion correction terms are a priori not hermitian. Furthermore, they are total energy 
dependent, like the unperturbed potentials. As explained in the introduction, these two difficulties are 
not independent. 



6 Conclusions 

We have written a 3D equation for three fermions (our basic three-cluster equation) by combining the three 
two-body potentials obtained by an exact 3D reduction of the corresponding two-fermion Bethe-Salpeter 
equations, putting the spectator fermion on the mass shell. In this way, the cluster-separated limits 
are still exact, and the Lorentz invariance / cluster separability requirement is automatically satisfied, 
provided no supplementary approximation, like the Born approximation, is made. This equation can 
be written in terms of the two-body potentials, or in terms of two-body transition operators (Faddeev 
formalism). The use of positive free-energy projectors in the chosen reductions of the two-fermion Bcthc- 
Salpeter equations prevents our 3D three-fermion equation from continuum dissolution. The potentials 
are hermitian and depend only slowly on the total three-fermion energy. The one high-mass limits of 
our "basic three-cluster equation" are approximately exact. The correction of the remaining discrepancy 
would demand the introduction of higher-order three-body terms. 

Our combination of cluster separability and Lorentz invariance in the three fermion problem makes 
explicit use of the fact that the clusters can only be two-fermion and/or free fermion states. This is not 
directly adaptable to four and more fermion systems. In this respect, 3 is still not N. 

The difference between our equation and Gross' spectator model equations consists in higher-order 
three-body contributions. Supplementary investigations would be necessary to decide which approach 
provides the best approximation to the three-fermion Bethe-Salpeter equation. 
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A lot of work remains to be done in order to precise a possible perturbation calculation program, 
correcting the remaining discrepancies with the three-fermion Bethe-Salpeter equation. We succeeded in 
deriving our 3D equation from an approximation of the three-fermion Bethe-Salpeter equation, in which 
the three-body kernel is neglected while the two-body kernels are approached by positive-energy instan- 
taneous expressions, with the spectator fermion on the mass shell. The correction terms are thus known 
at the Bethe-Salpeter level and can be transformed into corrections to the 3D equation, unfortunately 
not hermitian a priori. 

There exists an infinity of ways of performing the 3D reduction of the two-fermion Bethe-Salpeter 
equation. The potentials generated by these reductions could all by used to build a three-fermion 3D 
equation, keeping however in mind the continuum dissolution problem. This leaves us a large freedom to 
suit phenomenological needs. 

Our two-body potentials are the sum of an infinity of contributions symbolized by Feynman graphs. 
Keeping only the first one (Born approximation) or a finite number of them renders the Lorentz covariance 
of the two-fermion clusters only approximate. A Born approximation preserving the Lorentz invariance / 
cluster separability property can be obtained by using another 3D reduction based on a covariant second- 
order two-body propagator of Sazdjian, combined with a covariant substitute of A ++ . This leads to a 3D 
three-cluster equation which is covariantly Born approximable, but more complicated |]37j| . 
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